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Abstract 

In the last years it has been revealed that if measuring relative to higher 
order event planes higher order flow coefficients Vn for n > 2 can be 
measured. It also turned out that Bose-Einstein (HBT) correlation radii also 
show 3rd order oscillations if measured versus the third order event plane 
^ 3 - In this paper we investigate how these observables can be described via 
analytic hydro solutions and hydro parameterizations. We also investigate 
the time evolution of asymmetry coefficients and the mixing of velocity field 
asymmetries and density asymmetries. 


1 Introduction 

In relativistic heavy ion collisions, an expanding and cooling medium is created, 
usually referred to as the strongly interacting quark gluon plasma. Hydrodynamics 
provides a tool to investigate the time evolution of this medium, and exact analytic 
models are particularly useful in this regards. Usually spherical, axial or ellipsoidal 
symmetry is assumed in these solutions, as these are simple to handle and rep¬ 
resent geometries that yield realistic results for many soft observables. However, 
event-by-event fluctuating nuclear distributions yield event-by-event fluctuating ini¬ 
tial conditions, and thus higher order azimuthal asymmetries arise. In order to 
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describe these one has to utilize higher order asymmetries in hydro as well. This 
was successfully done in numerical calculations, see for example mm- 

In this paper we discuss the first exact analytic solutions [3] of relativistic hydro¬ 
dynamics that assume higher order asymmetries and thus give realistic higher order 
flow coefficients. We also discuss possible extensions of this approach, by analyzing 
the time evolution of the asymmetries in a numerical approach, and by investigating 
the effect of spatial versus momentum space anisotropies. 


2 Multipole solutions and higher order anisotropies 


The first 1+3D relativistic solution with ellipsoidal geometry was discovered in 
Ref. [4]. In this solution the thermodynamic quantities at a given proper-time r are 
constant on the surfaces of an expanding ellipsoid, defined by the s scale variable 


s = 


X- I _ ^ ^ 

' y2 ' ^2 ’ 


( 1 ) 


where rx,ry,rz are the spatial coordinates, while X, Y, Z are the time-dependent 
axes of the ellipsoid. The velocity profile as a function of space-time coordinates 
is given in form of a 3D Hubble flow, i.e. /r. With these, u^d^s = 0 holds 

(if the expansion of the axes is linear in time). In Ref. [4j it was already indicated, 
that more complicated scale variables can also be written up (with u^d^s = 0 
still holding). In Ref. [3] we showed that this solution can indeed be extended to 
multipole symmetries with a generalized scale variable 


rr.N 

s = ;^(l + ecos(A^(/))) (2) 

With the s given in Eq. the new solutions can be given in cylindrical coor¬ 
dinates (r, 0, z) as: 


® H cos(A^((/) - tAjv))) + ^ 

(3) 

n = nf v{s) 

(4) 


(5) 

/ Tf\ 3+3/K 

P = Pf[-) 

(6) 


and still representing a Hubble-flow, as in the original paper of Ref. [4]. In 
the formula for s, being the Nth order reaction planes (which cancel from the 
observables). This way we get new solutions with almost arbitrary shaped initial 
distributions, see Fig. It is important to note here that however, the initial 
state fluctuation in the observed collision is present through the orientation of the 
Nth order reaction planes and the strength of higher order asymmetries, the event 


2 



62=0.5, 83=0, 84=0 


82=0.5, 83=0.3, 84=0 


82=0.5, 83=0.3, 84=0.2 



Figure 1: Heat map of s values in the transverse plane, with multiple superim¬ 
posed symmetries. The more components are included, the more asymmetric 
the shape gets. 


plane orientation itself does not affect the measured quantities. Thus if every vn is 
measured relatively to the A^th order reaction plane, then the (event-through-event) 
averaged value of will correspond to an average n-pole anisotropy 

We also calculated hadronic observables from the above solution (see details 
of the freeze-out scenario in Ref p] or Ref. [5]. A comparison to PHENIX data 
on higher order harmonics measured in 200 GeV Au+Au collisions |6] is shown in 
Fig.H Fit parameters of the model are cn (for N = 2,3,4), ut and b (To and tq 
was fixed to values given from spectra and HBT comparisons of a similar model, as 
described in Refs. mm- 

3 Time evolution of the anisotropies 

In the above described solution, the anisotropies don’t change over time - due to 
the lack of pressure gradients and the Hubble-flow. In a numerical framework, we 
investigated how the introduction of pressure gradients, various speeds of sound 
and viscosity coefficients influence the time evolution of the asymmetries, when 
starting from an initial condition that is very similar to one described by know 
analytic solutions - except in pressure, where we used a pressure profile similar 
to the density profile given in usual Hubble-expansion models mm- We used a 
multi-stage predictor-corrector method outlined in Ref. [7]. This is a finite volume 
scheme, where the initial flux is a weighted average of the Lax-Friedrichs and Lax- 
Wendroff fluxes, called GFORCE. This flux is used to make a new prediction on 
the grid points, which is used to get a better flux approximation. This procedure is 
repeated for a number of times, as described e.g. in Ref. [7]. This multi-stage flux 
gives results that are comparable to those of the Godunov method. We tested our 
method with known analytic solutions given in Refs. PHE]. We analyzed both 
non-relativistic and relativistic hydrodynamics, and arrived at similar conclusions. 

Fig. shows the result of a nonrelativistic calculation of the time evolution of 
the energy density. If we assume a small amount of viscosity, it makes the flow 
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Figure 2: Fits to PHENIX 200 GeV Au+Au data [6] in 5 centrality bins. Fit 
parameters are summarized in Ref. [3] 
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t=Ofm/c t=0.5 fm/c t=0.76fm/c t=lfm/c 



t=Ofm/c t=0.5 fm/c t=0.76fm/c t=l fm/c 



Figure 3: Time evolution of energy density. Top row shows viscosity free case, 
the bottom row with /i = lOMeV • fm viscosity. 


Energy density spatial asymmetry coefficients 



Figure 4: The time evolution of the asymmetry coefficients 61 ^ 2 , 3,4 in the energy 
density, as modified by a small amount of viscosity (left plot) and the change in 
speed of sound (right plot). 
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Figure 5: The time evolution of the asymmetry coefficients in the energy density 
in a relativistic calculation, as modified by the change in k = c~^ (right plot). 


itself and thus the disappearance of asymmetries slower. The time evolution of the 
asymmetries themselves is shown in Fig. In this figure we also see the effect of 
speed of sound in this nonrelativistic calculation: the reduction of speed of sound 
makes the asymmetries disappear slower - due to the reduction of the speed of 
sound waves. A similar effect is seen in case of a relativistic calculation, as shown 
in Fig. the increase of k = c~‘^ makes the disappearance of asymmetries slower. 
It is important to see that this also slows down the speed of cooling, which means 
the system will freeze out latter. 

4 Anisotropy mixing 

It is important to see that there may be asymmetries in both momentum space 
(i.e. in the velocity field) and in density (i.e. in energy density or pressure), and 
both influence the measured anisotropies. To investigate this effect, we created a 
multipole version of the Buda-Lund model with a scale variable given in Eq.[^ 
but we introduced a multipole flow field as well. We start from a “flow potential” 
which gives us the flow: 


v = (7) 

2 

The flow field at a given time is spherically symmetric if ^ ^ with H being 

a Flubble-coefficient at that given time. Elliptical symmetry is obtained with $ = 
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^(1 + X 2 cos{2cp)), while 

2 3 

^ ^ (1 + X 2 cos{2ip)) + ^X3 cos(3(^)) (8) 

represents a triangular perturbation of the elliptical flow. Of course the various 
anisotropies can have various event planes (symmetry planes), but the specific angle 
of these does not enter into the results. Inspired by Ref. [10], we analyzed how X 2,3 
and 62,3 influence flow coefficients V 2 and vs - see results in Fig. Compared to 
the spatial anisotropy, velocity field anisotropy has a much larger effect on elliptic 
and triangular flow coefficients. 



X2 


X3 
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Figure 6: Dependence of flow anisotropy coefficients ^2,3 on asymmetry param¬ 
eters X2,3 and 62,3. Velocity field asymmetry has a stronger influence on flow 
coefficients. 


5 Summary and acknowledgments 

In this paper we showed an extension of the scope of analytic relativistic hydro¬ 
dynamics to higher order azimuthal asymmetries, compatible with realistic (event- 
by-event fluctuating) geometries. Higher order flow observables were calculated 
from this model, and are found to be compatible with data. In the analytic model, 
anisotropy parameters were independent of time, thus we investigated their time 
evolution in a numerical framework, developed for this purpose. We also inves¬ 
tigated how velocity- and density-field anisotropies “mix”, in the framework of a 
multipole Buda-Lund model. We are thankful to Tamas Csorgo and Marton Nagy 
for useful discussions with respect to this project. We thank the WPCF commu¬ 
nity and the WPCF 2014 organizers, in particular the local hosts, Tamas Novak 
and Tamas Csorgo, for the possibility to present this work. We also thankfully 
acknowledge the support of the OTKA grant NK 101438. 
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